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ABSTRACT 

We investigate the nature of high-z host galaxies of long Gamma- Ray Bursts (LGRBs) 
by means of state-of-the-art numerical simulations of cosmic structure formation and 
evolution of galaxies. We combine results from different runs with various box sizes and 
resolutions. By assigning to each simulated galaxy the probability to host a LGRB, 
assumed to be proportional to the mass of young stars, we provide a full description of 
the physical properties of high-z LGRB host galaxy population. We find that LGRBs 
at z > 6 are hosted in galaxies with typical star formation rates SFR ~ 0.03 — 
0.3 M yr _1 , stellar masses M* ~ 10 6 - 10 s M Q , and metallicities Z ~ 0.01 - 0.1 Z Q . 
Furthermore, the ratio between their doubling time and the corresponding cosmic 
time seems to be universally equal to ~ 0.1 — 0.3, independently from the redshift. 
The distribution of their UV luminosity places LGRB hosts in the faint-end of the 
galaxy luminosity function, well below the current capabilities of space- or ground- 
based optical facilities. This is in line with recent reports of non-detection of LGRB 
hosts using extremely deep HST and VLT observations. In conclusion, high-z LGRBs 
are found to trace the position of those faint galaxies that are thought to be the major 
actors in the rc-ionization of the Universe. 
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1 INTRODUCTION 

Long Gamma-Ray Bursts (LGRBs) are powerful flashes of 
7-rays that are observed with a frequency of about one 
per day over the whole sky. The 7-ray emission is accom- 
panied by a long-lasting tail, called afterglow, usually de- 
tected over the whole electromagnetic spectrum. Their ex- 
treme brightness easily over-shines the luminosity of their 
host galaxy and makes them detectable up to extremely 
high redshi fts, as shown by the discovery of GRB 090423 
at z = 8.2 (|Salvaterra et al.ll2009l: iTanvir et alj|2009h . or of 
GRB 090429B at z ~ 9.4 jCucchiara et al.ll201ll ). Metal ab- 
sorption lines can often be identified in their afterglow spec- 
tra, allowing a study of the properties of the environment 
in which they blow. Once the afterglow has faded, follow-up 
searches of the LGRB host galaxy become possible. 

Observations of high-z LGRBs can prov ide unique infor- 
matio n about the high-redshift Universe (see lMcQuinn et al.l 
120091 . and references therein). For example, LGRBs can be 
used to measure the neutral-hydrogen fraction of the inter- 
galactic medium (e.g. [McQuinn et al.ll2008l ; iGallerani et al.l 
120081 : iNagamine et al.ll2008h. to study the galaxy metal and 
dust content (e.g. lSavaglio et al.ll2012l ; IWang et alll2012i). to 
probe the intergalactic radiation field (e.g. Inoue et al.ll201ol ) 



or even to constrai n the level of primordial non-Gaussianities 
(Maio et all |2012| ). The tight observational link between 
high-z distant galaxies and LGRBs is particularly evident 
when looking at their conditions during the first few bil- 
lion years of life of the Universe. In fact, for the first time, 
surprisingly high metallicities (with super-solar abundances) 
have been estimated from GRB 090323 at redshift z ~ 3.5, 
in correspondence of a merging event between two galaxies 
jSavaglio et al.ll2012l ). The significant star formation rate 
(~ 6M0/yr) enhanced by their interaction might probably 
have triggered the birth of the LGRB progenitor and the 
simultaneous determination of a disturbed, metal-rich sys- 
tem, identifiable as a high-redshift, massive, sub-millimeter 
galaxy. 

Moreover, LGRBs might represent the most promising way 
to directly detect the very first stars: the so-called Popula- 
tion I I I (Pop III) star s (e.g.lCiardi fc Loebll20o"ol;lToma et all 



l201ll : ICampisi et ail l2011al : Ide Souza et alj l201ll ). These 



form in pristine environments and a re predicted to have 
large masse s (~ 1 - 10 2 M R ) (e.g IWooslev et ail |2002| : 
iHeger et all 120031 : IWooslev fc Hegerl 2006), because the 
most relevant coolants in such environments (H-derived 
molecules, like H2 or HD) are not able to cool the gas 
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efficiently and induce small-scale fragmentation. This is, 
instead, possible only in metal-enriched gas , where com- 
mon Population II-I (Pop II-I) stars can form l|Baraffe et al.l 
2001), either via lo w-temperature fine-structure transition 
meta l cooling (e.g. ISchneider et alj |2002| ; iBromm fc Loebl 
2003, when a minimum critical metallicity Z cr n ~ 10 4 Z© is 
reached), or via dust cooling (with a minimu m critical metal- 
licity as low as Z cr u ~ 1O _6 Z0, according to lSchneider et all 
12001 120121 ). The different initial mass functions (IMF) for 
Pop III stars can influence the thermodynamic state of the 
cosmic medium and the evolution of the following Pop II 
regime. Indeed, larger-mass stellar populations have harder 
spectral energy distribut ion (SEP) and can instantaneously 
produce more photons l|Schaererll2002l ). efficiently heating 
the cosmic gas. Also, they die and enrich the surrounding 
medium earlier, leading more rapidly (in any case, within 
~ 10 7 — 10 s yr) to the following Pop II star formation regime 
(e.g. iMaio et al-llioToh . Furthermore, since different progen- 
itor stars (i.e., massive Pop II or P op III stars) can poss i- 
bly affect the resulting LGRB rate (|Campisi et al.ll2011al ). 
to investigate the LGRBs properties it is crucial to have a 
good description of early structure formation, and to prop- 
erly model the transition from the pri mordial, Pop HI regime 
to the following, s tandard. Pop II one ijTprnatore et al 1120071 : 
IMaio et alJlioiol . l201ll ; IWise et alJl2012T ). 

In this paper, we will use state-o f-the-art N-body, hy- 
drodynamic, chemistry simulations l)Maio et al.l 120101 ) of 
structure formation and evolution of high- z galaxies to un- 
veil the nature and the physical properties of Pop II LGRB 
hosts at early epochs. Our work extends toward higher red- 
shifts the results already obtained by previous theo retical 
studies of the nature of LGRB hosts (ICourtv et al 



Nuza et all 120071: ICampisi et al.l I2OO9I : IChisari et al 



2004 



2010 



Campisi et al.1 2011bf ) that have mainly focused on the ex- 



pected properties of these galaxies at low or intermediate 
redshift where models can be compare directly with observa- 
tions. In spite of the different assumptions and approaches, 
all these works predict that LGRB should be hosted in low- 
mass, young, sta r-forming galaxies, c onsistently with obser- 
vations (see e.g. ISavaglio et al.l 120091 ). As this galaxy pop- 
ulation is believed to dominate the star formation activity 
in the early Universe |Salvaterra et al. Il201ll ; iJaacks et al.l 
l2012l ; lMunoz fc Loebll201ll ; lLacev et alj|201ll ), we would like 
to investigate to which extend LGRB can be used as a new 
tool to study the first phases of structure formation. 

Furthermore, extensive searches of high-z LGRB 
host galaxies has been rece ntly carried out with HST 
i|Tanvir et alj|2012h and VLT (|Basa et alj|2012h . In spite of 
the very deep limits reached by these observations, none of 
the targets have been identified. The lower limits on the UV 
luminosities reached using the deepest observations avail- 
able correspond to magnitudes Muv ~ —17 -I — 20, as pre- 
cisely reported in Table [1] together with upper limits for the 
star formation rate (SFR). These limits suggest that high- z 
LGRB hosts have lower SFRs with respect to their low-z 
counterparts |Basa et al.lfeOld '). By comparing available ob- 
servation limits to the results of our numerical simulations, 
we will study the detectability of high- z LGRB hosts with 
current and future instrumen ts. We will include in our anal- 
ysis GRB 060927 at z = 5.47 l|Ruiz-Velasco et alj|2007h . and 
the four known LG RBs at z > 6: GRB 0509 04 at z = 6.3 
l|Kawai et alj|200d ). GRB 080913 at z = 6.7 (|Greiner et all 



GRB 


z 




SFR [M Q /yr] 


060927 


5.47 


> -18.02 


< 0.65 


050904 


6.3 


> -19.95 


< 4.1 


080913 


6.7 


> -19.00 


< 1.3 


090423 


8.2 


> -16.95 


< 0.38 


090429B 


~ 9.4 


> -19.65 


< 2.4 



Table 1. High- 2 LGRB host observations from Tanvir et al. 
(2012) which are considered in this work. The columns refer to, 
from left to right: name of the LGRB, its redshift z, UV magni- 
tude yijjv and star formation rate. 



2009h. GRB 0904 23 at z = 8.2 l|Salvaterra et al.1 120091 ; 
Tanvir et all 120091 ) . and GRB 090429B w ith a photometric 
redshift of z ~ 9.4 jCucchiara et al.|[201lh - see also Table [fl 
We will also provide the expected LGRB host distributions 
as function of stellar masses (M*), gas phase metallicities 
(Z), and specific SFR (sSFR). This will give us a complete 
theoretical picture of the typical environments where Pop II 
LGRBs should be found at redshift z ~ 6 — 10, as inferred 
from the analyses of our numerical simulations. Regarding 
LGRBs powered by massive Pop II I stars, we refe r the in - 
terested reader to our previous work. [c"ampisi et al.l (|2011al ). 
where we compute their rate and the properties of their host 
galaxies by means of the same numerical simulations used 
here. 

The paper is organized as follows: in Section [5] we 
present the N-body, hydrodynamical, chemistry simulations 
used in this work; in Section [3] we compare the simulated 
host properties with available data; and, finally, we summa- 
rize our results in Section [4] Throughout the paper mag- 
nitudes are given in the AB system, logarithms are al- 
ways base-10, and the underlying cosmological model as- 
sumed for all the calculations and simulations is a stan- 
dard ACDM model, with total-matter density parameter 
at the present fio,m = 0.3, cosmological constant density 
parameter ilo.A = 0.7, baryonic-matter density parame- 
ter Qo,b — 0.04, expansion rate in units of 100 km/s/Mpc 
h — 0.7, spectral normalization as = 0.9, and primordial 
spectral index n = 1. Unless differently specified, metallici- 
ties refer to the gas phase. 



2 SIMULATIONS AND ANALYSIS 

In the next sections, we will briefly describe the essential 
features of our simulations (Section 12.11) . and then we will 
present the data analysis (Section l2.2l) leading to the results 
shown in the following Section [3] 



2.1 Simulations 



Here we use the numerical simulations bv lMaio et all (|2010h . 
in which gas collapse and condensation through atomic or 



molecular cooling (e . g. iGalli fc Palla j_998j : Ab e J e t al.l 



2002; 



2011 



Yoshida et all I200& IMaio et al l |200d . 120071 , 120091 . 
Wise fc Abel 120071 : IWise et al.l l2012h is followed down to 
proto-galac t ic scales, where star fo r mation activity (e.g 



Katj Il992l; iHerncmist fc Katzl [l989l : iHerncmist fc Springe! 
20031 : ISpringel fc Hernqu ist 2003) takes place and induces 
feedback effects on the ambient medium. Supernova (SN) 
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Model 


L [Mpc/ft] 


rrigas [M Q //i] 


m DM [M©/ft] 


r\ [kpc/ft] 


MSFRIMeyr- 1 ]) 


L30 


30 


9 x 10 6 


6 x 10 7 


4.7 


> 


L10 


10 


3 x 10 5 


2 x 10 6 


1.0 


[-1-5,0] 


L5 


5 


4 x 10 4 


3 x 10 5 


0.5 


< -1.5 



Table 2. Simulation parameters. The columns refer to, from left to right: model name; comoving box length, L; gas particle mass, m, 
dark matter particle mass, r«DM; softening length, rj\ star formation rate, M*. 



explosions dist ribute metals in the surrounding gas (as de- 
scrib e d in e .g. Springel fc Hernquistl 120031 ; iTornatore et al.l 
|200o1 . I2007I . l20ld ) altering its chemical com position on the 
short timescales typical of massive stars (|Wooslev et al.l 
l2002i ; lHeger fc Wooslevll2002l ; ISchaererll2002l ). 

More in detail, for each stellar partic le we follow 



JWooslev fc Weaverl [l99l 
Thielemann et alj [200^ 



its timescales (Padovani & Matteucci 199 and yields 

Heeer fc Wooslevl 12001 l20ld: 



van den Hoek fc Groenewegenl 



19971 ) according to a given IMF (see later). When SNe 
take place the surrounding gas is randomly 'kicked' with 
a velocity of ~ 500km/s (kinetic feedback) and heated up 
at temperatures of roug hly ~ 10 5 K (thermal feedback). 
Simultaneously, different metal species (C, O, Si, Fe, Mg) 
are spread around the neighbouring regions according to the 
corresponding yields, by mean of a kernel-based smoothing 
procedure that mimics metal diffusiorQ. Dependencies on 
numerical resolution of the differ ent sub-grid mo d els ar e 
minor, as a lso st udied in, e.g., iMaio fc Iannuzzl (|201ll ); 
iMaio et alj |201ll ). and they usually affect only the very 
first bursts of star formation, to converge relatively soon 
at slightly later times. Further effects of different IMFs, 
SN rages, yields determinations, or possible changes in the 
b asic SPH implemen tation have been discussed extensively 
in lMaio et ail |2010l ). 

The transition between the Pop III and Pop II star for- 
mation regime mentioned in the Introduction is assumed to 
happen at a critical metallicity Z cr u = 10 _4 Zq, such that 
if the star forming gas has Z < Z cr u, then a Salpeter-like 
Pop III IMF is assumed with mass range between 100 Mq 
and 500 Mq a nd a pair- instability S N range between 140 Mq 
and 260 M Q l|Wooslev et all [2002); otherwise, a standard 
Salpeter IMF is adopted with mass range between 0.1 Mq 
and 100 M , and SNII range between 8 Mq and 40 M . 
The cosmological field is sampled for dark matter and bary- 
onic matter at redshift z — 100, following a common Gaus- 
sian matter distributiorQ . In particular, here we consider 
three cubic volumes with comoving sides of 30 Mpc/h, 
10 Mpc//i, and 5 Mpc/h. The identification of the simulated 
objects (with their gaseous, dark and stellar components) 
was carried out by applying a friends-of-friends technique, 
with comoving linking length of 20% the mean inter-particle 
separation, and a minimum number of 32 particles. To en- 
sure and resolve properly the low-mass end of our simulated 



We note that metal diffusion is a complex issue and currently 
cannot be properly treated neither in SPH codes (because of the 
difficulty in diffusing the material above the smoothing scale), nor 
in grid codes (that are affected by over- mixing problems). 
2 Deviations from Gaussia nity are not expected to significantly 
alter our conclusions (see IMaio fc Iarmuzzil 1 20 111; iMaiol l201ll ; 
IMaio fc Khochfar|[20l3) . 



objects, we will consider only well-resolved structures with 
at least 100 baryon particleqj (the minimum number of dark 
matter particles per object is usually an order of magnitude 
larger). Moreover, we check that our results do not change 
significantly if only objects with more than 300 baryon par- 
ticles are considered. Substructures are identified by using 
a SubFind algorithm, which follows a spherical-overdensity 
approach and additionall y discriminates am ong bound and 
non-bound particles (e.g. iDolag et al.| [2009). The parame- 
ters used in the simulations are listed in Table d] For more 
detai ls we refer the reader to the original paper (|Maio et alj 
l2010h . 

Before discussing the analysis of the simulations, it should 
be noticed that, despite the amount of detailed physi- 
cal implementations included, some ingredients are still 
missing or could be further improved. For example, metal 
diffusion is accounted for by smoothing over the SPH 
kernel: this approach might work reasonably on small 
scales, but it will probably have some limitations on 
large scales. Neither radiative feedback (e.g lGnedin fc Abell 
200ll;IWhalen fc Normanll2008h nor thermal conduction (e.g. 



Dolag et al.ll2004l ) are included, and both are likely to intro- 



duce some modifications in the properties of the host galax- 
ies. On the other hand, other specific issues, like numerical- 
viscosity schemes or smoothing kernels to refine the SPH 
algorithms are likely to give major improvements when re- 
solving the details of hydrodynamical shocks and turbulence, 
but they s hould not a ffect significantly the general trends we 
find here (|Pricdl2012l . and references therein). 

2.2 Simulation analysis 

Rare galaxies with very high SFRs are properly accounted 
for only in the larger simulation box, while galaxies at very 
low SFR are resolved only in the small simulation box which 
has a better numerical resolution. Therefore, we jointly con- 
sider the results of different simulation boxes by selecting 
galaxies on the basis of their SFRs, which are reported in 
the last column of Table [2] The SFR thresholds are iden- 
tified by inspecting the number density distributions and 
correspond to the SFR bins in which the galaxy number den- 
sity is consistent in the larger and smaller box. We check at 
which level both numerical resolution and the unavoidable 
lack of small er, unresolved, star fo rming structures in the 
larger boxes (jNagamine et alj 120041 ) alter our results going 
from the smaller to the larger simulation box. We find that, 



3 I Bate fc Burker t (1997) have shown that objects containing 3-4 
times the number of neighbors on which density is computed are 
properly resolved. In our simulations, 32 neighbors are considered 
so that the cut at > 100 baryon particles is a reasonable choice 
to avoid numerical problems. 
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Figure 1. Luminosity functions at z = 6, 7, 8 and 9. Long-dashed 
(green), short-dashed (red) and dotted (blue) histograms show the 
result of L5, L10 and L30, respectively, while solid (black) his- 
tograms refer to the luminosity functions obta ined using jointly 
the three simul ations. Data ar e from : iBouwens et al.l J2007I. cir- 
cles) a t 2 = 6; IBouwens et al] ||20H|. ci r cles) and lMcLure et ah! 
i2010l . triangles') at z = 7. 8; lQesch et all ||2012| . squares) at z = 8. 
The cont i nuos line at z = 6 is the fit to the LF as obtained by 
ISu et al] fcOllT) . 



in the mass-SFR plane, SFRs roughly agree when consid- 
ering objects with the same stellar masses extracted from 
the three different boxes. However, we obtain a clear, con- 
stant displacement in the gas phase metallicities of galaxies 
with the same mass (or SFR) going from L 5 to L10 and L30 
sugge stive of a numerical resolution effect ()Nagamine et al.l 
2004). Therefore, we renormalize the metallicities of simu- 
lated objects by adding 0.1 dex and 0.5 dex to all galaxies 
(at all redshifts) in L10 and L30, respectively. We notice 
that since most of the LGRBs (~ 90%) form in objects se- 
lected from L5 and L10, the rather large correction applied 
to metallicities in L30 will not affect significantly our final 
results (see later). 

The total luminosity of a galaxy at wavelength A, L\ 
(in erg s" 1 Hz -1 ), forming Pop II stars at a rate SFR 11 and 
Pop III stars at a rate SFR 111 , is then given bjQ 



t u ,Z) SFR n + /i n SFR ln r m , 



(1) 



where 



l x (r , Z) is the SED template (in 
erg s -1 Hz -1 Mg 1 yr) corresponding to the energy 
emitted per unit time and frequency (erg s -1 Hz -1 ), and 
per unit star formation rate (in Mq yr -1 ) for Pop II stars 
with mean age r 11 and metallicity Z. Similarly, I™ is the 



We notice that, despite their similar notations, in the following 
and H n are dimensionally two different quantities. 



SED template (in erg s _1 Hz -1 Mg 1 ) corresponding to the 
energy emitted per unit time and frequency (erg s _1 Hz -1 ), 
and per unit stellar mass (in Mq) fo r Pop HI stars with 
mean lifetime r IH = 2.5 x 10 6 yr jSchaererl I2OO2T ). We 
implicitly assume that the emission properties of Pop 
III stars are roughly constant during their short lifetime. 
Finally, we neglect dust abs o rption in our calculation. 
As shown by ISalvaterra et al.l |201ll ). dust extinction for 
galaxies with M t < 10 s M at z > 6 is E(B — V) < 0.01, 
while a larger contributi on is expected in m ore massive and 
metal-rich objects (e.g. iDaval et al] l20ld ). As it will be 
clear in the next Section, the large majority of LGRBs are 
hosted in low-mass galaxies so that our assumption of no 
dust will not affect our conclusions. 

As a consistency check, in Fig. [T] we show the lumi- 
nosity function (LF) of galaxies at z = 6, 7, 8 and 9. 
The solid histogram reports the LF as obtained by con- 
sidering jointly the three simulation boxes. The resulting 
LF is compared with data from the literature. In the range 
of luminosities in which observational data and simulation 
results overlap, we find a good agreement with the am- 
plitude of the observed galaxy LF at all redshifts. More- 
over, our_sjmulations,simik numeri- 
cal ((Salvaterra ct al. 2011; Jaacks et al. 2 0121) a nd semi- 
analytical (M ufioz fc Loebll201ll ; Lacev et al] 2013 ) studies, 
predict a rather steep faint-end slope consistent wit h the 
estimates based on HUDF data (jBouwens et al.ll201ll ). The 
agreement is even more remarkable as no attempts have been 
made to fit or adjust the theoretical curves to the observed 
LF, i.e. they have been computed directly from the simula- 
tion output with no free extra parameters. We note that the 
bright-end of the LF at z = 6 is slightly overestimated by 
our simulations, in particular by L30. This effect is indeed 
expected since, as mentioned above, we neglect dust extinc- 
tion important for bright massive objects. We note however 
that this does not affect the conclusions of our paper, as 
less than 5% of LGRBs are found to be hosted in galaxies 
brighter than Muv ~ —20 (see Fig. [3]). In conclusion, our 
simulations provide a good description of the galaxy popula- 
tion in the early Universe in the luminosity range probed by 
HST and match the extrapolation of the observed faint-end 
slope in the redshift range we are interested in. In the follow- 
ing, we will make use of these results to predict the nature 
and the properties of LGRB host galaxies in the redshift 
range z = 6 — 10. 



3 SIMULATED LGRB HOST GALAXIES 

Long GRBs are believed to be li n ked to the death of mas - 
sive stars (e.g. iFrver et al. 1 1 19991 : IWooslev fc Hegerl l2006t ). 
This idea is supported by the association of local and low-z 
LGRBs with the explosion of a broad-line type Ic SN (e.g . 
Delia Valle et al]|2003l; iHiorth et al"1l2003l ; iPian et al]|2006l ; 
Bufano et alJl2012l ; lMelandri et alj|2012h . While the precise 
physical details of the LGRB phenomena are still matter 
of debate, it is reasonable to assume that the probability 
for a LGRB to be hosted in a given galaxy should be pro- 
portional to the mass of young stellar particles (with ages 
below ~ 10 7 yr, corresponding to stellar masses > 30 Mq) 
in each galaxy. Here we neglect the LGRBs powered by Pop 
III stars, whose contribution is likely to be sub-dominant 
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Figure 2. Differential probability of finding LGRBs per logarithmic SFR bin in L30 (left panel), L10 (central) and L5 (right). Each 
color pixel refers to a galaxy and is color coded according to the probability in eq. l2t . 



with respec t to the Pop II popu lation in the redshift range 
considered l|Campisi et aljfioilal ) . This choice is further mo- 
tivated by the fact that observed high-z LGRBs are quali- 
tatively similar to lo w- and intermediate-z o nes suggesting 
common progenitors (|Salvaterra et all I2009T ) . Additionally, 
in selecting galaxy hosts we do not c onsider any metallicity 
bias for the LGRB formation. Indeed. lCampisi et al.l |2011bl ) 
have shown that strong cut-off in the progenitor metallicity 
of the order of 0.1 — 0.3 Zq are at odd with the consistency 
between LGRB host galaxies at z < 1 and the observed Fun - 
damental Metallicity Relation (FMR; lMannucci et aTfeoilh . 
While higher critical metalli cities for the LGRB progenitors 
(as su ggested for example bv lNuza et al.ll2007l ; lGeorgv et al.l 
120091 ) can not be excluded at the present stage, we want 
to stress that our results are not affected by the existence 
of such critical threshold. Indeed, the large majority of the 
young stellar particles from which we compute the probabil- 
ity of the galaxy to host a LGRB event have stellar metallic- 
ities Z < 0.3 Zq in the redshift range considered. Therefore, 
our results do not represent in any way a test for the exis- 
tence of a metallicity bias for the LGRB formation. 



where, d-/VG.RB(log(SFR)) is the number of LGRB expected 
in each SFR bin (in units of M0/yr), normalized by their 
total number, Ngrb- 

We display the resulting distributions in Fig. [2] for L30 
(left panel), L10 (central) and L5 (right). In general, the 
SFRs vary over a wide range of values (see discussion in 
Sec. l3.3p . Typically, in L30 more massive objects with larger 
star formation rates are found, while in L5 small objects are 
resolved which are not present in the other two boxes. Inter- 
mediate values of the SFR are found in L10. From the Fig- 
ure, it is clear that the most frequent values for the LGRB 
probability defined in eq. ([2]) are around or below ~ 0.15 
(blue) in L30, compared to ~ 0.15 — 0.30 (between blue and 
yellow) in L10, and > 0.30 (yellow and light yellow) in L5. 
This shows that, contrary to what naively expected, LGRBs 
are more probably found in intermediate- or small-size ob- 
jects with lower SFRs, and not in (often highly star forming) 
large systems, due to their paucity in the cosmological land- 
scape. This is particularly relevant for high- redshift detec- 
tions and investigations of primordial stars and also suggests 
that LGRB host luminosities should be quite dim (see next 
Section). 



3.1 LGRB host distribution and SFR 

A quantitative estimate of the probability of finding a LGRB 
in a given galaxy can be simply obtained studying the 
SFRs. At a first glance, one would think that larger objects 
have larger SFRs and hence a larger probability of form- 
ing LGRBs. However, while this holds for a single galaxy, 
this statement is not necessarily true on a global scale, as it 
does not take into account the lower cosmic abundance of 
bigger structures. Thus, we measure the statistical probabil- 
ity of LGRB hosts by considering a sample formed by the 
three simulations described in Section [2] and computing the 
differential probability distribution of LGRBs, APgrb, per 
logarithmic SFR: 

APgrb djV G fl B (log(SFR)) . . 

dlog(SFR) Ngrb dlog(SFR) ' 1 ' 



3.2 Luminosities of high-z LGRB host galaxies 

The distribution of UV luminosities of our simulated hosts 
is shown in Fig. [3] together with the l imits obtaine d in th e 
deepest HST/WFC3 observations by iTanvir et al.l (|2012l ). 
As explained in previous Sections, the sample is obtained by 
jointly considering the results of our three simulation boxes. 
At z = 6 the probability shows a wide distribution, but 
only a very small fraction of LGRBs would reside in galaxies 
brighter than Muv = —18, i.e. accessible to current instru- 
ments. This reflects the fact that most of the young stars 
form in faint galaxies. At higher redshifts, the peak of the 
distribution shifts towards fainter galaxies, and mimics the 
evolution of the characteristic luminosity of the galaxy LF. 
In particular, at z — 8 only ~ 10% of LGRB hosts would be 
brighter than the extremely deep limit obtained for LGRB 
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Figure 3. Distributions of LGRB host absolute magnitude Mjjy 
at z = 6, 7, 8, and 9, as indicated by the legends. Arrows show 
the lower limits obtained by Tanvir et al. (2012), with the cor- 
responding LGRB name. In each panel, the solid line shows the 
differential distribution (scale on the left y-axis), while the dotted 
line shows the normalized cumulative distribution (scale on the 
right y-axis). 



090423. The non-detection of LGRB hosts at z > 6 is then 
supporting the idea that most of the star formation arises in 
the faint-end of the galaxy LF at those early epochs, when 
mini-haloes host the initial phases of primordial star forma- 
tion. Our results also support the hypothesis that LGRB 
host galaxies are likely to be too faint to be detected by 
currently available space- and ground-based facilities. Deep 
JWST observations are therefore needed to image high-z 
LGRB host galaxies. We find that in order to have a 50% 
probability to detect a LGRB host at z — 6 at 5a level, an 
observation of ~ 2 x 10 5 sec should be foreseen with the 
F150W filter on NIRCAIvfl For z = 8 LGRBs, ~ 3 times 
larger integration times are expected. 

These results suggest that LGRBs can be used as sign- 
post for those very faint, but extremely common, galaxies 
that are now thought to reionize the Universe. This is further 
illustrated in Fig. [4j where we show the fraction of ionizing 
photons produced in galaxies fainter than a given absolute 
magnitude M,vv- Most of the ionizing photons are provided 
by galaxies much fainter than those resolved by HST/WFC3 
observations and will be resolved only in the deepest JWST 
fields. The limits on the luminosities of high- z LGRB hosts 
are over-plotted in the Figure. Less than 10% of the ioniz- 
ing photons are provided by galaxies brighter than the lim- 
its reached by current LGRB host searches. This holds true 
even in the case of GRB 090423, in spite of the extreme limit 
reached by the dedicated HST observations. Thus, high-z 
LGRBs provide the first (best) targets for JWST and ELT 



5 On the base of current version of the JWST exposure time 
calculator, http://jwstetc.stsci.edu/etc/ 




Figure 4. Fraction of ionizing photons produced in galaxies 
fainter than Mjjy. Arrows show the lower limits obtained by 
Tanvir et al. (2012) with the corresponding GRB name. Solid 
(dashed) line refers to z = 6 (z = 8) galaxies. 



observat ions aimed to seek the sources of cosmic reionization 
(see also lTrenti et alj|2012h . 



3.3 Statistical properties of high-z LGRB host 
galaxies 

Now, we further explore the other physical properties of 
LGRB host galaxies in the redshift range z = 6 — 10, for 
our whole simulated sample. In Fig. [5] from top to bottom, 
the panels show the distribution of SFRs, stellar masses, 
specific SFR (sSFR, i.e. the SFR per units of stellar mass) 
and gas phase metallicities. In all panels the dotted lines 
refer to the cumulative fraction. 

The distribution of SFRs indicates that most of the sim- 
ulated hosts have SFR ~ 0.03 - 0.3 M yr" 1 with only a 
few percent having SFR > 1 Mq yr . This again is in line 
with the limits derived from indivi dual observat i ons o f the 
fields of the most distant LGRBs bv lTanvir et alj (|2012h and 
iBasa et al l l|2012i ). Indeed, considering current observational 
limits, the probability of non detection for all the five targets 
is as high as 70% when the SFR distribution of simulated 
galaxies is adopted. It is interesting to compare this value 
to the probability of non detection expected assuming that 
2 > 6 galaxies have a distribution similar to that of LGRB 
hosts detected at z < 1. Basa et al. (2012) report that in this 
case the probability is only ~ 6%, suggesting that the SFR 
distribution of high-z hosts should peak at much lower values 
than low-z one, as indeed found in our simulations. Assum- 
ing that LGR Bs are good tracers of the formation of stars 
(but see e.g. ISalvaterra fc Chincarinill2007l ; ISalvaterra et all 
2012), current observational limits may already provide the 
evidence that at high z the bulk of stars resides in objects 
with much lower SFRs than in the local Universe, as also ex- 
pected from numerical simulations of early structure forma- 
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Figure 5. Properties of simulated LGRB hosts in the rcdshift 
range z = 6 — 10. Panels from top to bottom show the distribu- 
tion of SFR (in M yr" 1 ), stellar masses (in Mq), specific SFR 
(in Gyr ), and gas phase metallicities (in solar units, Z/Zq). 
The dotted line shows the normalized cumulative distribution (see 
right y-axis). The histogram with arrows in the top panel reports 
the limits on the SFR of LGRB hosts as obtained by Tanvir et al. 
(2012). The data point in the bottom panel refers to the mctal- 
licity inferred from Si absorption lines in the afterglow spectrum 
of GRB 050904 at z = 6.3, by Kawai et al. (2006). 



tion (e.g. I Abel et alJliool lYoshida et ai1l2003l : | Maio et al. l 
l20ld ; IWise et all 12012ft . 

The majority of the simulated LGRB hosts is found to 
have relatively low stellar masses, M* ~ 10 6 — 10 8 Mq. As 
a result, in spite of their relatively low SFR, high- 2 LGRB 
hosts have a very large sSFR, of ~ 3 — lOGyr -1 , showing 
that they are experiencing a strong burst of star formation. 
The mean sSFR of high- z hosts is much larger than the 
observed low- 2 LGR B one (estimated to be ~ 0.8 Gyr -1 ; 
ISavaglio et al.l [2009V However, their doubling time (equal 



to the inverse of the sSFR), tdb, is ~ 0.1 — 0.3 of the Hubble 
time, tn{z), similar to what found for the low- 2 LGRB hosts, 



tdb(z) _ 1 

t H (z) ~~ sSFR(z) t H (z) 



0.1-0.3 



(3) 



at all redshifts. 

In spite of the high z, LGRB hosts have gas phase metal- 
licities in the range —0.5 < log(Z/Z ) < —3, with a peak 
in the distribution around Z ~ 0.03 Zq. These values are 
in line with those inferred from absorption line measure- 
ment^] in the LGRB afterglow at lower redshift, 2 < z < 4, 
suggesting that the properties of the environment in which 
LGRBs explode are not so different in the early Universe. 
This scenario is supported by the metallicity inferred from 
the analysis of the afterglow spectrum of GRB 050904 at 
z = 6.3 (as marked in the bottom panel of Figure (5jQ. 
Besides that, no other metallicity constraint at z > 6 has 
been obtained so far owing to the low signal-to-noise ra- 
tio of the spectra. However, we note that other LGRBs de - 
tected at z = 4-5, e.g. GRB 090205 ilD'Avanzo et al.l l 2010t ) 
and GRB 100219A (|Thoene et all 12012ft . show metallicities 
Z ^ 0.1 Zq, further supporting our findings. This is also 
in line with numerical studies of cosmic metal enrichment, 
that predict efficient pollution from early s tellar populations 
llTornatore et al.1 120071 ; iMaio et all |2010| . l201ll ; Wise et all 
|2012j), and the consequent presen ce of significative am ounts 
of metals in high-z LGRB hosts (jCampisi et al.|[2011al ). 



3.4 Physical relations for LGRB host galaxies 

To further characterize the environment properties of high- 
z LGRBs, we compute the luminosity-metallicity relation 
(M[/v-Z), mass-metallicity relation (M-Z) and FMR rela- 
tion of the LGRB hosts. These are shown in Fig. [6] for 
z — 6 and z = 8, where contour levels refer to the loci 
in which we expect to find 30%, 60%, and 90% of the sim- 
ulated LGRBs. We find that LGRB hosts, and simulated 
high-z galaxies in general, follow well defined relations in 
these planes. The luminosity-metallicity relation shows the 
largest scatter that reduces significantly considering the M- 
Z relation. An even tighter relation is found when the pa- 
rameter fj,o.32 = log(M*) — 0.32 log(SFR) is considered. This 
trend is similar to what foun d for field galaxies at lower red- 
shift ilMannucci et all 120101 ) and for LGRB hosts at 2 < 1 



IjMannucci et al.l 120111 ). As expected, LGRBs are tracing 
the low-mass faint-end of the relations. Measurements of 
metallicities in high-2 LGRB afterglows can therefore pro- 
vide unique information about the metal enrichment history 
of very small galaxies in the primordial Universe. 

Adopting the Muv-Z relation and considering the ob- 
servational limits on the LGRB host luminosities, we can 
infer that the bursts detected so far at z > 6 should have 
exploded in galaxies with Z < 0.1 Zq. Moreover, we find 



6 Metal absorption lines over-imposed on the LGRB afterglow 
probe properly the metal content along the line-of-sight and not 
necessarily of the galaxy itself 

7 A slightly lower metallic i ty for GRB 050904 has been computed 
recently by iThoene et al.l 1120121 ), log(Z/Z Q ) = -1.6 ± 0.1. This 
value is still consistent with our findings. 
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Figure 6. Luminosity-mctallicity (left columns), mass-metallicity (central columns) and Fundamental Metallicity Relation (right 
columns) for the LGRB host galaxies at z = 6 (upper rows) and z = 8 (lower rows). Contour plots report the 30%, 60%, and 90% 
probability of hosting a LGRB. Arrows refer to Tanvir et al. (2012) results and, in the absence of a measured metallicity, have been 
positioned arbitrarily at Z = 0.3 Zq, while the metallicity of LGRB 050904 has been obtained by Kawai et al. (2006). 



that the host galaxy of GRB 050904, in order to be con- 
sistent with the simulated relation, should be at least two 
magnitudes fainter than current upper limits. Thus, its de- 
tection would result very difficult even in extremely long 
HST observations as done in the case of GRB 090423. 

iThoene et alj l|2012T ) presented a sample of LGRBs at 
2 = 4 — 5 for which the measure of the metallicity and of the 
host galaxy luminosity has been obtained. F our out of the six 
LGRBs considered by IThoene et al.1 ((2012) show metallici- 
ties > 0.1 Zq. The corresponding LGRB host luminosities in 
the rest-frame B band are in the range [-22,-20]. The other 
two LGRBs have Z < 0.1 Zq, but only upper limits on the 
host luminosity are given, with Mb ^ —21. Although these 
bursts lie at redshifts slightly lower than those studied here, 
we note that the position of these bursts in the L-Z plane is 
consistent with our findings. In particular, they support the 
idea that LGRB host galaxies should be already enriched at 
the level of a few percent of the solar metallicity at z 6. 

It is interesting to compare these maps with the analo- 
gous fo r Pop III LGRB hosts, as computed bv lCampisi et all 
feoilaT ). With respect to Pop III LGRB hosts, we find that 
Pop II LGRB hosts are more massive and more metal en- 
riched at all redshifts. Indeed, most of the Pop III LGRBs ex- 
plode in galaxies near the critical metallicity line. Although 
some is found to be hosted in the outer edge of more evolved 



galaxies - see Figure 3 in lCarrrpisi et al.1 (|2011al ). Low metal 
abundances and peculiar abundance ratios may help in dis- 
criminating between Pop III and Pop II LGRBs. 



4 CONCLUSIONS 

We investigated the nature of high-z host galaxies of LGRBs 
by means of N-bo dy, hydrodynamical, chemistry simulations 
l|Maio et al.ll2010T ) of galaxy formation and evolution, trac- 
ing both Population III and Population II-I star forming 
regimes, and providing a good description of the general 
field galaxy population at z > 6 in the luminosity range 
probed by current HST observations. 

We predicted the physical properties of the LGRB host 
population in the light of our understanding of galaxy forma- 
tion and evolution in the early Universe. We found that high- 
redshift (z ~ 6 — 10) LGRBs are hosted in bursty galaxies 
with typical star formation rates SFR ~ 0.03 — 0.3 Mq yr _1 , 
stellar masses M* ~ 10 6 — 10 8 Mq, specific star formation 
rates sSFR ~ 3 — lOGyr -1 (larger than the low-z counter- 
parts of up to one order of magnitude), and metallicities 
Z ~ 0.01 — 0.1 Zq. Furthermore, the ratio between the re- 
sulting LGRB host doubling time and the corresponding cos- 
mic time seems to be universally equal to roughly ~ 0.1 — 0.3, 
independently from the redshift. 
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The distribution of their UV luminosity places simu- 
lated LGRB hosts in the faint-end of the galaxy luminos- 
ity function, well below the current capabilities of space- or 
ground-based optical facilities, in line with recent reports of 
non-detection of LGRB hosts in extremely deep HST and 
VLT observations. Deep JWST observations should be fore- 
seen in order to detect these objects. These results suggest 
that LGRBs are the signposts of those faint, but extremely 
common, galaxies that are thought to rcionize the Universe. 
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